Photon number resolution using time-multiplexed single-photon detectors 
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Photon number-resolving detectors are needed for a variety of applications including linear-optics 
quantum computing. Here we describe the use of time-multiplexing techniques that allows ordi- 
nary single photon detectors, such as silicon avalanche photodiodes, to be used as photon number- 
resolving detectors. The ability of such a detector to correctly measure the number of photons for 
an incident number state is analyzed. The predicted results for an incident coherent state are found 
to be in good agreement with the results of a proof-of-principle experimental demonstration. 
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I. INTRODUCTION 

There has been considerable interest recently in the 
development of photon detectors that are capable of re- 
solving the number n of photons present in an incident 
pulse. Photon number-resolving detectors of this kind 
are needed for a linear optics approach to quantum com- 
puting 0,0) for example, and they may have other appli- 
cations as well 0, such as conditional state preparation 
0, H . Here we describe a simple time-multiplexing tech- 
nique that allows ordinary single-photon detectors to be 
used as photon number-resolving detectors. The achiev- 
able performance of these devices is comparable to that 
of cryogenic devices 0, 13 that are being developed 
specifically for number-resolving applications. Detectors 
based on an atomic vapor have also been proposed [9tll0|. 

Photon number-resolving detectors can be character- 
ized in part by the probability P(m\n) that m photons 
will be detected in a pulse that actually contains n pho- 
tons. For linear optics quantum computing, it is essen- 
tial that P(n\n) be as close to unity as possible. As 
one might expect, this requires the single-photon detec- 
tion efficiency r\ to be as large as possible in order to 
maximize the probability of detecting all of the photons 
present in a pulse. Silicon avalanche photodiodes (APDs) 
operating in the visible have a relatively large value of r\ 
and a small dark count rate, but they only produce a 
single output pulse regardless of the number of incident 
photons. The time-multiplexing technique described here 
avoids that limitation while retaining the large value of 
r\ and other potential advantages associated with the use 
of commercially- available APDs. 

The basic idea of the time-multiplexing technique is 
to divide the incident pulse into N separate pulses of 
approximately equal amplitude that are displaced in time 
by a time interval At that is larger than the detector dead 
time r. The probability that one of the divided pulses 
contains two or more photons becomes negligibly small 



for N ^> n, in which case the fact that an APD can only 
produce a single output pulse is no longer a limitation. 

Our implementation of such a time-multiplexed detec- 
tor is illustrated in Figure ^ An incident pulse propa- 
gating in a single-mode optical fiber was divided equally 
into two paths using a 50/50 coupler. The difference L 
in path lengths was sufficiently large that the propaga- 
tion times differed by more than the detector dead time. 
The two pulses were then recombined at a second 50/50 
coupler and split equally into two more paths differing 
in length by 2L, which results in four pulses all sepa- 
rated in time by At. This process was repeated once 
more with a path-length difference of 4L, after which a 
fourth 50/50 coupler split the pulses again and directed 
them into one of two silicon APDs. The net result was 
the creation of TV = 16 pulses of approximately equal 
amplitude. Unlike earlier approaches involving detector 
arrays 0, 0, 0, 0] , this technique requires only two 
detectors, and it avoids the need for optical switches and 
a storage loop [IE EE ■ 

The remainder of this paper is organized as follows: 
In Section [nj we present a theoretical analysis of the 
expected performance of a time-multiplexed detector 
(TMD) as a function of the detector efficiency and pho- 
ton losses. The analysis makes clear that high detector 




* michacl.fitch@jhuapl.edu 
1 bryan.jacobs@jhuapl.edu 
' todd.pittman@jhuapl.edu 
-james.franson@jhuapl.edu 



FIG. 1: Implementation of a time-multiplexed detector 
(TMD) for photon number-resolving applications. An inci- 
dent pulse is divided into N weaker pulses separated by a 
time interval At using a series of fiber couplers and optical 
fiber delay lines. For sufficiently large values of N and At, this 
allows two ordinary silicon avalanche photodiodes to measure 
the number of incident photons with a probability of success 
that depends on their detection efficiency. 
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efficiency is needed for any photon number resolving de- 
tector, including the TMD and various cryogenic detec- 
tors under development jEOlOl- Section ITTT1 describes an 
experimental implementation of a TMD and the results 
obtained using coherent state inputs. A summary and 
conclusions are presented in Section llVI 



II. THEORETICAL ANALYSIS 

In the most likely application of a time-multiplexed de- 
tector, an incident pulse will contain an unknown number 
n of photons and the TMD is intended to measure the 
value of n. Since the measurement process destroys any 
coherence between states with different values of n, it 
is sufficient to consider an incident number state (Fock 
state) with a specific value of n and then calculate the 
response of the detector in the form of the probability 
distribution P(m\n). The response to an incident coher- 
ent state will be considered in the next section. 

It is therefore assumed that a number state \n) = 
(l/Vn!) (a t )"|0) is incident on the initial 50/50 coupler, 
where the operator a' creates a photon in the incident 
mode and |0) is the vacuum state. The effects of the first 
coupler can be described by the operator transformation 

at-(ot + aJ)/V2 (i) 

where the operators a\ and a, create photons in the 
shorter and longer paths, respectively. Similar trans- 
forms can be applied to represent the remainder of the 
50/50 couplers. After these transformations, the state 
vector will involve n photons distributed over 16 differ- 
ent modes. 

It is important to include the effects of photon loss, 
which will degrade the performance of the device in a way 
that is similar to the effects of limited detector efficiency. 
The largest source of loss is absorption and scattering in 
the optical fiber. This can be described by the fraction 
/ of the incident power that is transmitted through a 
length L of fiber. (In the experiments described in the 
next section, / has an approximate value of 0.97). The 
transmission through the fiber section of length 2L is 
then / 2 , etc. In order to simplify the analysis, it was 
assumed that the short sections of fiber have negligible 
length and that losses in the fiber connections are also 
negligible. 

Losses in the optical fiber were included in the analysis 
by inserting an additional fiber coupler into each of the 
longer loops. The reflection and transmission coefficients 
of these couplers were adjusted to give a transmission 
probability of / (or the appropriate power of /) for each 
of the sections of fiber. An additional field mode corre- 
sponding to the photons removed from the loop by these 
couplers for each of the pulses was included in the state 
vector, which increased the total number of modes from 
16 to 23. The effects of photon loss could then be taken 
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FIG. 2: The calculated probability distribution P(m\n) for 
observing m detection events given that n photons were in- 
cident on a time-multiplexed detector. These results corre- 
spond to the case of two incident photons (n = 2) and a 
detector efficiency r] = 0.7 with no loss (/ = 1). 



into account by applying operator transformations anal- 
ogous to Eq. ^ 

After all of the operator transformations had been ap- 
plied, the state vector contained a large number of terms 
corresponding to all of the ways in which n photons can 
be distributed over the 23 modes. As a result, mathe- 
MATICA was used to sum the contribution of each term in 
the state vector to the probability distribution P(m\n). 
In doing so, it was assumed that Po = (1 — r/) 9 , where Po 
is the probability that a detector will detect no photons if 
q photons are incident upon it. (This ignores any possi- 
ble correlations between the effects of multiple photons.) 
The probability Pa that a detector will detect at least 
one photon and produce an output pulse is then given by 
Pa = 1 — Po- Having determined the values of Po and 
Pa for each of the detectors, it was straightforward to 
calculate the probability of obtaining m detection events 
and add that contribution to the probability distribution 
P(m\n) for each term in the state vector. 

This analysis includes the fact that each of the 16 
pulses will have different probabilities of reaching the de- 
tector, since they all travel through different lengths of 
optical fiber. In the limit of no loss (/ — > 1) and perfect 
detection efficiency (77 — > 1), the TMD is mathematically 
equivalent to a lossless multiport device, whose output 
properties can be obtained analytically [12j . Our numer- 
ical calculations agree with the analytic results in that 
limit. 

Figure |21 shows the results of this analysis for the case 
of two incident pho tons (n = 2), which is a case of inter- 
est pi 0, E, IT^.l20| in linear optics quantum computing. 
The detection efficiency was taken to be rj = 0.7, which is 
typical of commercially-available silicon APDs, while the 
losses were assumed to be zero (/ = 1) in this example. 
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FIG. 3: Calculated probability distribution for observing m 
detection events when two photons are incident and the de- 
tection efficiency is 0.2 with no loss (/ = 1), as in Fig. [5] 



The limited value of N = 16 reduces the probability of 
detecting the correct number of photons to 15/16, while 
the limited detection efficiency further reduces the prob- 
ability of a successful measurement to ~45%. It can be 
seen that a detection efficiency of 0.7 substantially limits 
the ability of a TMD to resolve the number of photons 
present in a pulse, even for the simple case of n = 2, and 
higher detection efficiencies will be required in order to 
improve the performance. 

Similar results are shown in Figure |31 for the case of 
r) = 0.2, which is comparable to the estimated external 
efficiency of the superconducting detector in Ref. @ . De- 
tectors of that kind are mathematically equivalent to tak- 
ing the limit of large N, which would increase the prob- 
ability of obtaining the correct photon number by 6% 
(1/16) above that shown in Figure|3| Nevertheless, it is 
apparent that further increases in the detection efficiency 
will be required. 

Larger numbers of incident photons give a correspond- 
ingly smaller probability of measuring the correct value 
of n. As an example, Figure 01 shows the probability dis- 
tribution P(m\n) for the case of n = 5, / = 0.97, and 
r\ = 0.43. These parameters correspond to the exper- 
imental apparatus described in the next section, where 
the reduced value of rj includes the combined averaged 
effects of input and output coupling loss, connector and 
splice losses, and excess losses in the fiber couplers. 

More generally, the probability P(n\n) of measuring 
the correct number of photons is plotted in Figure [S] as 
a function of the detector efficiency. For simplicity, the 
losses were assumed to be zero here, in which case P(n\ri) 
has the analytic result 

16' 

P(n n) = — j— rr Tj n for n < 16 (2) 

v 1 ' 16 n (16-n)! ' - \ j 

as shown in the Appendix. In addition to the effects of 
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FIG. 4: Calculated probability distribution for m detection 
events for a Fock state with n — 5. Here the parameter values 
were r\ = 0.43 and / = 0.97 which correspond to the experi- 
mental conditions of Section HTT1 when the additional losses in 
the connections and couplers are included. 
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FIG. 5: Calculated probability P(n\n) of detecting all n inci- 
dent photons as a function of the effective detection efficiency 
r\. The points were calculated numerically using the operator 
transformation technique of Eq. while the lines correspond 
to the analytic formula of Eq. [5] 



limited detector efficiency, it can be seen that the limited 
number of secondary pulses (N — 16) also has a signif- 
icant effect as the value of n increases, even for perfect 
detector efficiency. At telecom wavelengths, the losses in 
optical fibers are much smaller and it should be possible 
to use much larger values of N. 
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III. EXPERIMENTAL RESULTS 



In order to perform a proof-of-principle demonstra- 
tion of a time-multiplexed detector, we implemented the 
TMD of Figurenusing single- mode optical fiber and 2x2 
fiber couplers designed to give 50/50 reflection and trans- 
mission at a wavelength of 702 nm. Since a source of 
number states is not currently available for large values 
of n, the experiments were performed using a weak co- 
herent state. This was intended to allow a comparison of 
the theory with experiment and a demonstration of the 
basic feasibility of the TMD approach. 

The input and output ports of the TMD had FC-style 
connectors while all of the other connections were made 
using a fusion splicer in order to minimize the losses. 
The two detectors were commercial silicon APD single- 
photon detectors (Perkin-Elmer model SPCM-AQR-I3) 
with measured deadtimes of approximately 60 ns. The 
length L of the first long loop was chosen to be 22 m, 
which gave a delay between adjacent pulses of At — 
110 ns, which was substantially larger than the dead- 
time. Short pulses (50 ps duration) were generated us- 
ing an externally triggered fiber-coupled diode laser at 
680 nm. A custom-made electronic circuit was used to 
count the total number of photon detections in a speci- 
fied time window following each trigger. The triggering 
rate was kept sufficiently small that there was no overlap 
of the sequence of N pulses from two different triggering 
events. 

The optical fiber used for the delay loops had a 4 /im 
core, which supported single-mode propagation of visible 
light (610 nm to 730 nm). The attenuation at 700 nm 
was specified by the manufacturer to be 6 dB/km, which 
corresponds to a transmission factor of / = 0.97. Small 
but significant losses also occurred in the input and out- 
put couplers, fusion splices, and excess loss in the 2x2 
fiber couplers (typically 3% each). The total transmis- 
sion of the fiber multiplexing device was measured to be 
0.55, averaged over all N pulses. 

The photon number distribution of the coherent state 
at one of the outputs of the TMD is given by a Poisson 
distribution: 



P(li,n) 



(3) 



where \i is the mean photon number. For a perfect de- 
tector {r\ = 1), P(fi,0) would give the probability P of 
detecting no photons. It can be shown from the linear- 
ity of the system that this result can be generalized to 
the case of an imperfect detector by replacing \i with 
fi' = t)\l. (The limited detection efficiency is equivalent 
to placing an attenuator in front of the detector.) Thus 
Po = cxp(-fi') and Pa = 1 — Pa, where Pa is defined 
once again as the probability that a detector will detect 
at least one photon and produce an output pulse. 

For values of / close to 1 , the losses for all N pulses can 
be taken to be approximately the same without signifi- 
cantly affecting the distribution of the total number of 
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FIG. 6: Histogram of the number of events in which m pho- 
tons were detected for an incident coherent state pulse with 
a relatively large intensity. The light bars correspond to the 
theoretical prediction of Eq. [I] based on a least-squares fit 
with rjljio = 13.1, while the dark bars correspond the the 
experimental measurements. 



counts. In that case, // = rjl[i /l% for each of the pulses 
incident on the detectors, where the transmission factor 
I accounts for the average total loss in the system (mea- 
sured value I = 0.55), rj is the APD detection efficiency 
(typ. 0.7), and /i Q is the mean number of photons inci- 
dent upon the TMD. The probability V(m) of obtaining 
exactly to detection events (for m < 16) is then given by 
the binomial distribution 



V(m) 



16! 



(16 — to)! to! 



,16— m 



(l-Po) 



(4) 



A comparison of the theoretical predictions of Eq. 0] 
with the experimental results from a relatively intense 
coherent state are shown in Figure The data corre- 
sponds to a histogram of the total number of events in 
which to photons were counted in an arbitrary data col- 
lection interval, which is proportional to the probability 
distribution V(m). The theoretical results correspond to 
a least-squares fit to Eq. 21 with two free parameters, a 
normalization constant (which reflects the length of the 
data collection interval) and the value of fjf, which cor- 
responded to a best-fit value rjlfi = 13.1. It can be seen 
that the data are in relatively good agreement with the 
theoretical predictions. The small discrepancies between 
the theory and experiment are probably due to the ap- 
proximation that all 16 pulses are subject to the same 
loss. 

Similar comparisons of the theoretical predictions and 
the experimental results are shown in Figure \7\ for a less 
intense coherent state (r?i/x = 2.65) and in Figure[S]for a 
relatively weak pulse (rjlfio = 0.57). It can be seen that a 
TMD can provide reliable estimates of the mean number 
of photons in a coherent-state pulse, but it can also be 
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FIG. 7: Comparison of the theoretical prediction and exper- 
imental results for an incident coherent state with a smaller 
intensity than in Figure [§] Here the best fit corresponds to 
r)lfj, = 2.65. 
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FIG. 8: Comparison of the theoretical prediction and experi- 
mental results for an weak coherent state input, as in Figure 
[(J Here the best fit corresponds to r)l[/, = 0.57. 



seen that the observed peak in the data is substantially 
less than the true mean number of photons. 



IV. SUMMARY 

In summary, we have described the use of time- 
multiplexing techniques to allow ordinary silicon 
avalanche photodiodes to be used as photon-number de- 
tectors. Although silicon APDs have high efficiencies 
and low dark counts, they only produce a single out- 
put pulse regardless of the number of photons that are 
incident. This difficulty can be avoided by splitting the 
incident pulse into a time sequence of N weaker pulses, 



provided that N 3> n. The use of time multiplexing does 
lengthen the effective detector response time, which may 
also be of importance depending on the intended appli- 
cation eiei m m. 

A proof-of-principle experiment was performed using 
optical fiber loops for the case of N — 16 and a coherent 
state input pulse. Good agreement was observed between 
the theoretical predictions and the experimental results. 

The calculated response of a TMD to an incident num- 
ber state makes it clear that high detection efficiencies 
will be required in order to truly resolve the number of 
photons in an incident pulse. But the same comment ap- 
plies to any other detector intended for photon number- 
resolving applications, such as the cryogenic detectors 
that are being developed for that purpose. The results 
of our analysis suggest that the most important require- 
ment for true number-resolving detectors is a high single- 
photon detection efficiency; an intrinsic ability to resolve 
photon numbers may be of lesser importance, since that 
capability can also be achieved using ordinary single- 
photon detectors. 

Finally, we would like to note that a similar device 
has been independently demonstrated and analyzed by 
Achilles et al. 24]. 

APPENDIX: 

In this Appendix, we provide a derivation of Eq.[2J and 
generalize it to arbitrary numbers of detected photons. 

Non-classical photon statistics have been observed [2{| 
at a multiport beam splitter when the path lengths are 
matched, and to describe such effects, a quantum formal- 
ism must be used (as in Section [H] above) . In contrast, 
the path lengths in the TMD are far from matched, with 
relative delays of > 100 ns. As this is two or three orders 
of magnitude longer than the coherence length, the pos- 
sibility of photon interferences can be neglected, and the 
relevant detection probabilities can be calculated from 
classical probability theory. 

Let the input mode be evenly distributed to N out- 
put modes, so that a single photon on the input has 
a probability 1/N of reaching any given output mode. 
We further assume there are no losses, corresponding to 
/ = 1 above, so this model describes a balanced, lossless 
Af-port. At each output mode, there is a detector with ef- 
ficiency rj. Suppose that n photons are sent to the input. 
Let the probability of m detection events (0 < m < n) be 
written P^(m\n). In the case of zero detections, P^ (0|n) 
is given by: 

if (Q|n) = (l-Tj) n (A.l) 

which assumes that each photon fails to be detected with 
an independent probability of 1 — r\. For the case where 
all photons are detected, P^ (n\n), must include a prob- 
ability of (rj/N) n for each of the ways in which a photon 
can be detected. Since all n photons must have gone 
to distinct output modes, there is a combinatoric factor 
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which counts the number of different ways to distribute 
n objects among TV bins, so that: 



P n N (n\n) = (| 



TV! 



(TV - n)\ 



for n < TV. (A.2) 



Eq. IA.2I reduces to Eq. for the special case in which 
TV = 16. 

For detectors of unit efficiency (rj = 1), Paul et al. [L2| 
gave a recursion relation for (in our notation) P^L 1 (m\n), 
along with a closed form solution. Their results can be 
extended to include non-unit detection efficiency by using 
the recursion relation: 



P?(m\n + l) = Pf(m\n) (1 - V ) + ^ 



TV 



+ P*(m-l\n) (N+l-m) 



TV 



(A.3) 



where the first term on the right hand side describes an 
additional photon failing to be detected or else exiting in 
a previously occupied output mode. The second term on 
the right hand side describes the detection of an addi- 
tional photon in a previously unoccupied output mode. 
In the limit r/ — > 1, the recursion relation Eq. IA.3I agrees 
with that given by Paul et al. (see Eq. 9 in Ref. 12]). 



Using the boundary conditions of Equations IA.1I and 
IA.2I the recursion relation Eq. IA.3I can be solved, with 
the result: 



P»(m\n) 



TV 



3=0 



(-iy 



(i-v) + 



(m - j)t] 



TV 



(A.4) 

for 77i < 7i < TV, where the binomial coefficient is ( N ) = 
TV!/[m!(TV -to)!]. The derivation of Eo. IA.4I is lengthy, 
but it can be verified that it satisfies the recursion rela- 
tion. 

Although Eqs. IA.1I through I A. 41 were derived using 
classical probability theory, they give exactly the same 
numerical results as the field operator approach of Sec- 
tion im 
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